IMEC21 0.001 AUS PATENT 
A METHOD AND APPARATUS FOR CHANNEL ESTIMATION 

5 This application claims priority from and incorporates by reference, in its entirety, 

U.S. Provisional Application No. 60/208,447, entitled "METHOD AND APPARATUS 
FOR CHANNEL ESTIMATION," filed May 31, 2001. 
FTRLD OF THE INVENTION 

The invention relates to modern telecommunication methods wherein equalization 

10 is exploited and channel estimation is needed. 

BACKGROUND OF THE INVENTION 
In telecommunication systems, wherein transmission from at least one 
transmitting peer to at least one receiving peer is considered, knowledge of the channel 
15 between said transmitting peer and said receiving peer, can be necessary for or be used 
for improving the determining of the transmitted signal at the receiving peer. Said signal 
improving process is also denoted equalization. Said channel knowledge can be obtained 
by performing a channel estimation, meaning determining the frequency response of said 
channel. 

20 In telecommunication methods based on subband processing one typically 

considers a channel with Nc carriers. Note that in subband processing of a data signal 
having a data rate, comprises in principle of splitting said data signal in a plurality of data 
signals, with a lower data rate and modulating each of said plurality of data signals with 
another carrier. Said carriers are preferably orthogonal. Channel estimation for channels 

25 with multiple carriers boils down to determining the channel frequency response at the 
frequencies of said carriers or at a subset thereof. Note that in modern wireless 
telecommunication methods the channel time domain response is short, thus having a 
finite time domain response, meaning a finite amount of non-zero samples. 

An example of a modern telecommunication method is Orthogonal Frequency 

30 Division Multiplexing (OFDM), which has become increasingly popular during the last 
decades, mainly because it provides a substantial reduction in equalization complexity 
compared to classical modulation techniques. Indeed, OFDM with cyclic prefix can be 
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equalized by a single low-rate complex multiplication on each carrier. For this reason, it 
has been adopted in the upcoming standards for high data-rate wireless networks, such as 
ETSI Hiperlanll and IEEE 802.11a. As opposed to former standards using OFDM 
modulation, the new standards rely on coherent QAM modulation and thus require 
5 channel estimation. Hence, the complexity of channel estimation is of crucial importance, 
especially for time varying channels, where this has to be performed periodically or even 
continuously. 

Channel estimation can be performed during reception of information or data 
between the transmitting and receiving peer, which is denoted blind channel estimation or 

10 by transmitting specific information, also denoted reference tones, between said peers. 
Although the invention focuses towards non-blind channel estimation, it is clear that a 
combination with blind channel estimation is possible. 

When relying on transmission of said specific information or reference tones, it is 
clear that it is beneficiary if one can estimate the channel frequency response at the 

15 frequencies of its carriers or at a subset thereof from only a limited set of reference tones, 
meaning less reference tones than the amount of carriers. In the context of multi-user 
transmission schemes, wherein each user transmits from a different transmitting peer 
towards at least one receiving peer, one is even obliged to use only a limited set of 
reference tones per user, if one wants to determining channel frequency responses of each 

20 of the users simultaneously. One can then for each user use a different set of reference 
tones, preferably scrambled all over the spectra of the channel. This scrambling is also 

denoted interleaving. 

Note that as non-blind channel estimation invokes transmission of specific 
information, which is in fact non-useful data from a users point of view, channel 

25 estimation is considered as an overhead, which should be efficient, meaning taking a 
small amount of time, thus a minimal amount of reference tones, but also an efficient 
implementation, meaning taking a small amount of operations in order to determine the 
channel estimate, resulting in a fast determination and taking a small area in chip 
implementations. However, although complexity reduction is aimed at, still optimality of 

30 the estimation operation is an important feature. 
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It must be emphasized that in the modem telecommunications standards, using 
multi-carriers approaches, the reference tones used for channel estimation can not be 
chosen to be equidistant as said standards define certain carriers to be zero carriers, 
meaning not to be used for transmission. Indeed the standards use some zero-carriers for 
5 spectral shaping, for example, some carriers are not used to allow smooth decaying of the 
spectral power on the border of the bandwidth. Hence, these standards (and possible 
variations) will be referred to as spectral shaping systems. 

Existing literature recognizes that, due to the structure of OFDM signals, the 
channel can be estimated by using the frequency correlation of the channel. This 
1 0 frequency correlation has inspired different approaches. 

Edfors et al. [O. Edfors, M. Sandell, J.J. van de Beek, S. K. Wilson, and P. O. Borjesson. 
"OFDM Channel Estimation by Singular Value Decomposition," IEEE Trans, on 
41 Communications}, 46(7):93 1-939, July 1998.] use explicitly the frequency correlation 

; 3 and derive a linear minimum mean squared (LMMSE) estimator. Using optimal rank 

;t 15 reduction, they develop a low complexity algorithm which computes an approximated 

LMMSE estimator. Complexity is here a trade-off with optimality. 
r Raleigh and Jones [G.G. Raleigh and V.K. Jones. ^Multivariate Modulation and Coding 

5 for Wireless Communication," IEEE Journ. on Special Areas in Communications, 

W 17(5):851-860, May 1999.] link the frequency correlation to the maximum delay spread 

S 20 and estimate the channel from a part of the carriers only. These carriers must be regularly 
*~ spaced, which limits the application of their method. 

It can be stated that prior-art channel estimation methods can not cope with non- 
equidistant reference or pilot tones, required by modern standard. Complexity reduction 
of the prior-art channel estimators leads to loss of optimality. 

25 

SUMMARY OF THE INVENTION 
It is a first aspect of the invention that a channel estimation method is presented 
with provides a maximum-likelihood estimate of the frequency response of the channel. 
A channel between at least one transmitting peer and at least one receiving peer is 
30 considered. Nu reference tones are transmitted. Said maximum-likelihood frequency 
domain estimation operation is such that the finiteness of the time domain response of the 
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channel is explicitly exploited, resulting in both an optimal but still a low complexity 
approach. The method can deal with multi-carrier channels with N c carriers and works 
also when the amount of reference tones Nu is less than N c . Note that with multi-carrier 
channel is meant a plurality of subchannels, each subchannel being centered around a 

5 carrier. Reference tones can be associated with a subchannel or a carrier. 

In an embodiment of the invention it is recognized explicitly that the determining 
of a maximum-likelihood estimate can be performed in the time domain. Therefore one 
recognizes a transformation of the received or captured reference tones from the 
frequency domain into the time domain, a time domain estimation step and a back 

1 0 transformation of the obtained time domain estimates to the frequency domain. 

In a further embodiment one recognizes that said transformations can be 
performed by applying IFFT and FFT operations, which can be efficiently implemented. 
Moreover the maximum-likelihood estimate determination is shown to be an efficient 
operation also, as it can be done by making linear combinations of the reference tones, 

1 5 requiring an amount of multiplications of the order of N H 2 , with N H the amount of non- 
zero valued time samples of the impulse response of the channel, which is typically small 
in subband processing telecommunication schemes such as OFDM used for wireless 
indoor communications. Note that said complexity reduction does not result in lack of 
optimality. The estimate of the channel retains to be a maximum-likelihood estimate. 

20 In an embodiment of the invention the applicability of the method also in case of 

non-equidistant reference tones is established. 

In an embodiment of the invention the exploitation of the fmiteness of the time 
domain response of the channel is made explicit by characterizing the dimensions of the 
operations to be performed, more in particular the subband and inverse subband 

25 processing, or in case of OFDM the FFT and IFFT operations. Also the weighting of the 
time domain measurements, being the transformed captured reference tones, used to 
obtain the maximum likelihood frequency domain estimates, is characterized by the 
fmiteness of the time domain response of the channel. 

In an embodiment of the invention the determining of the weighting operation via 

30 a weighting matrix is characterized as being as solving a least squares problem in the 
time domain in a non-iterative way, meaning that an explicit solution is presented. 
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In another embodiment of the invention the determining of the maximum 
likelihood estimate is characterized as a one step operation, which can be understood as 
solving a constrained least squares problem in the frequency domain, in a non-iterative 
way. In a further embodiment thereof, the method is characterized as a projection. 

In an alternative embodiment of the invention it is recognized that said method 
can be implemented as a two-step matrix multiplication process instead of a three-step 
approach. Indeed by splitting the weighting matrix in two submatrices, together defining 
said weighting matrix by multiplication, said submatrices preferably being equal, and by 
modifying the IFFT and FFT partial matrices, by left-hand and right-hand matrix 
multiplication with such a submatrix, the method can be executed as a two-step process 
wherein first an intermediate result is determined by multiplying a partial weighted N H by 
Nu IFFT matrix with said reference tones with N H the length of said finite time response 
and thereafter multiplying said intermediate result with a partial weighted N F by N H FFT 
matrix, to thereby obtain said maximum likelihood frequency domain estimates, is done. 

In a second aspect of the invention telecommunication circuits, including a 
channel estimation circuit, having the necessary subcircuits for executing one of the 
above described methods, are presented. 

In an embodiment of said second aspect a general purpose programmable circuit, 
programmed for executing said channel estimation methods, is disclosed. 

In another embodiment the channel estimation circuit, implements a three-step 
approach, each step having its dedicated circuitry. In a further embodiment thereof part of 
said dedicated circuitry comprises IFFT and FFT circuits. 

In another embodiment the channel estimation circuit, implements a two-step 
approach, each step having its dedicated circuitry. 

In another embodiment the channel estimation circuit, implements a one-step 
approach, each step having its dedicated circuitry. 

One aspect of the invention to proposes a channel estimation method, being used 
in both a single and multi-user context, which is capable of working with non-equidistant 
reference tones, and thus compatible with modern telecommunication standards. Further 
said channel estimation method should have low complexity without losing optimality of 
the channel estimate. 


BRTEF DESCRIPTION OF THE DRAWINGS 
Figure 1 relates to a modern telecommunication system exploiting a kind of 
subband processing (OFDM) wherein the invention can be used for performing 
5 equalization of the received signal. 

Figure 2 shows situation wherein the invented channel estimation outperforms 
classical approaches because it still provides an optimal estimate although no sufficient 
references tones (full line) are available to cover all carriers (dashed line) of the channel 
to be estimated. 

10 Figures 3-4 and 6-7 show comparison between the invention and other methods 

and between various embodiments of the invention. Comparisons of the FFT based and 
optimal methods of the invention with the SVD approximation methods (LMMSE) with 
respect to complexity and Bit-Error-Rate are shown. 

Figure 5 shows a typical topology as found in indoor wireless communication 
1 5 wherein the transmitting scheme with the invented estimation method can be exploited. 

Figure 8 shows an architecture of a signal transmitting and/or receiving circuit, 
wherein said channel estimation method fits. 

Figure 9 shows a three step approach digital circuit (10) for determining at said 
receiving peer from Nu received reference tones (100) said maximum likelihood 
20 frequency domain estimate (200) of said channel response at N F predetermined 
frequencies by directly exploiting the fmiteness of the time response of said channel. 

Figure 10 shows a digital circuit implementing a two-step approach comprises a 
first subcircuit (1020) transforming said received reference tones (100) into an 
intermediate results (400); and a second subcircuit (1040) transforming said intermediate 
25 result (400) into said maximum likelihood frequency domain estimates (200). 

Figure 11 shows a flowchart of the three step approach wherein reference tones 
(500), being frequency domain measurements, are used as input and a channel estimate 
(530) is obtained. The first step is a transformation from frequency domain to time 
domain (600), obtaining time domain measurements (510), the second step is an 
30 estimation step in the time domain (610), resulting in time domain estimates (620), 
thereafter being converted in a last step (620) into frequency domain estimates (530). 
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Figure 12 shows a flowchart of the two step approach wherein reference tones 
(500), being frequency domain measurements, are used as input and a channel estimate 
(530) is obtained. The first step (700) transforms said reference tones into an intermediate 
results (800) . The second step (720) further transforms said intermediate results into a 
channel estimate. 

Figure 13 shows a flowchart of a one step approach wherein reference tones (500) 
are transformed into a channel estimate (530) in a single step (900). 

Figure 14 shows a transmitting peer (1000) with at least one transmitting device 
(two shown) and a receiving peer (1010) with at least one receiving device or receiver 
(two shown). The channel estimation within said receiver tries to characterize the channel 
(1020) in between the transmitting device and the receiving device. 

DETAILED DESCRIPTION OF THE INVENTION 
The invention explicitly uses the finite delay spread of the channel and develops a 
low complexity algorithm capable of estimating the channel from part of the carriers 
only. A deterministic model is introduced and therefrom the associated Maximum 
Likelihood estimator is derived. This ML estimator can be interpreted as a transformation 
from frequency domain to time domain and back to frequency. The actual estimation is 
done in the time domain, where the number of parameters (i.e. the channel length) is 
small. The estimator is obtained by minimizing a quadratic criterion, which, combined 
with the small number of parameters, leads to a low complexity algorithm. As such, an 
exact low complexity solution is obtained. The invention is applicable when reference 
tones are available for all N c carriers of the multi-carrier channel but also when less 
reference tones are transmitted as in Pilot Symbol Assisted Modulation (PSAM). In an 
alternative formulation the invention can be casted as a Constrained Least Squares (CLS) 
solution. In such framework also explicitly the length of the channel impulse response is 
used and an ad-hoc constrained least squares estimator is derived. The invented method 
allows non regular spacing of the carriers. 

The invention is further described for OFDM but it is clear that any subband 
processing method can be dealt with. After introducing the OFDM system model, the 
statistical model and channel estimator is presented, along with the extension to PSAM 
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and the link with the CLS estimator. The time-frequency interpretation of the ML 
estimator is given, and the benefits provided by a combined PSAM/decision-feedback 
system are indicated. The complexity of the algorithm is presented and simulation results 
for indoor PSAM and spectral-shaping systems are discussed. 

5 Normal letters represent scalar quantities, boldface represent vectors and boldface 

capitals matrices. Slanted (resp. roman) letters indicate time (resp. frequency) domain 
quantities. X T , X H and X* respectively mean transpose, conjugate transpose and Moore- 
Penrose pseudo-inverse of X . If X is full column rank, then P x = X(X H X)" 1 X H is the 
orthogonal projection onto the space spanned by the columns of X . Further the invention 

10 refers to IFFT matrices (related to IDFT transformations) and FFT matrices (related to 
DFT transformation), which are defined as follows: 
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OFDM modulation consists in multiplexing QAM data symbols over a large number of 
20 orthogonal carriers. To this end, the QAM symbols of an OFDM symbol are passed 
through an Inverse Fast Fourier Transform (IFFT/TDFT). In the presence of a time 
dispersive channel, a Cyclic Prefix (CP) is prepended to each OFDM symbol to preserve 
orthogonality between carriers and eliminate InterSymbol Interference (ISI). A single 
user / single channel communication setup (see figure I), with OFDM modulation, is 
25 considered now, but it is clear that the estimation method can easily be extended to a 
multi-user, multi-channel context as known in spatial division multiple access systems. 
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Recall that the goal is to obtain h[0], h[N c -l], in order to perform channel 
equalization. The set -up is described by : 

Y = X©H+N 0) 


5 where the Hadamard (i.e. element-wise) product of the columns of X with H is used. For 
a single OFDM symbol this results in: 

y = x©HH-n 

The QAM source is written as X= [xo ... x M -i] , where x m = [x 0 , m ... x Nc .i, m ] is an 

10 OFDM symbol. N c denotes the number of carriers, m is a time index (often omitted for 
clarity) and M is the number of OFDM symbols. After IFFT and cyclic prefix insertion, 
the transmitted signal is X m = [ X Nc - N P ,m »• XNc-i.m Xo, m ... X N c-i,ni ] where N p is the 
size of the prefix and X= [ X 0 ... X M -il • For a channel h=[h 0 ... h Nh .i] , where N H being 
smaller than or equal to N p , the received vector is, after prefix removal and FFT, y m = x m 

15 H , where H= [H 0 ... Hnc-i] t is the FFT of the channel. Equations (I) and (2) further take 
the additive (possibly colored) Gaussian noise into account. Equalization is then done by 
a complex division on each carrier. 

Two types of training (PSAM and spectral shaping systems) are considered. With 
training here is meant determining or estimating of the channel frequency response. In 

20 classical training based estimation, all components of x are known. Spectral based 
systems use a minor modification of the classical training, zeroing a small number of 
carriers (named zero carriers) at the edges and in the middle of the band used. PSAM on 
the other hand bases it's channel estimation on a small fraction of the carriers, usually 
evenly spaced on the whole band, and possibly on varying positions from one OFDM 

25 symbol to the next. The invention is applicable in all said types of training. 

The ML estimator is derived, based on a reduced order model. It is extended to 
PSAM, linked to the CLS method and interpreted in terms of time-frequency 
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transformations. It is also applied on a combination of PSAM and decision-feedback 
estimation. 

First the case is explained wherein a channel with N c subchannels is determined 
by using an amount Nu of reference tones, also denoted pilot tones being equal to Nc.. 
5 The channel response is determined at the frequencies defined by its carriers or 
subchannels. 

Equation (2) shows that the OFDM system can be described as a set of parallel 
Gaussian channels. Since the time domain channel h has a finite length (in a well- 
designed OFDM system smaller than the prefix length) these parallel channels feature 
10 correlated attenuations. Considering, w.l.o.g., x= [1 1 ... 1] , equation (2) becomes 


y =F 

where F is a N c x N c FFT matrix, y is a Gaussian random variable with mean F [ h 
0] T and covariance C„ n . The signal y can also be denoted the captured reference tone at 
15 the receiving peer. The signal part of y is contained only in the space spanned by its 
mean. Separating the signal subspace from the noise only subspace, the received signal 
can be rewritten as : 

y = [F*F* 


h 
0 


+ n = H + n, 


(3) 


0 


4-n. 


(4) 


20 Relying on this, the reduced space signal is defined as 


r = Ffy = * +Fj[n = * +v, (5) 

where v is a zero mean Gaussian noise of covariance C vv = Fh H C nn Fj, . If C nn = o n 2 Inc , 
v is a white Gaussian noise of covariance matrix a a 2 l m • The reduced space signal has a 
25 log likelihood function expressed by 
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log/(r ) = -\ 0 tf%te(C w )) - (B*y- h fC^y-* ). (6) 


Maximizing this log likelihood with respect to h leads to the ML estimator given by 


(7) 


H = F fc Fty=P Ffc y ) 

5 where P F h denotes the orthogonal projection on the column-space of F h . Before 
performing the ML estimation, N H , the amount of non-zero time samples of the impulse 
response of the channel, must be determined, and we denote it's estimation as N H e- 

Formula (7) shows how the maximum likelihood frequency domain estimate of 
said channel response is obtained from said N H = N c captured reference tones y at said N F 

10 = N c predetermined frequencies. The same formula (7) shows that an orthogonal 
projection of said reference tones on the column space of a partial N c by N H FFT matrix 
is involved. 

Now the cases wherein less reference tones are used than subchannels in the 
channel is explained. In the case of Pilot Symbol Assisted Modulation and spectral 
15 shaping systems, not all symbols in x are known, and only a subset of N„ measured 
carriers or captured reference tones can be used. Only this part of the signal (denoted y u ) 
will be used. These are the reference tones. Equation (5) becomes 


r=^y tt = h^n = h+v, (8) 


20 where F has been decomposed as : 


F 



■ 




(9) 


and where measured pilots have been grouped together. Ungrouped pilots can be handled 
by straightforward permutations in columns and lines of the vectors and matrices. The 
25 ML estimator for spectral shaping systems, corresponding to (8) is H u = P Fu h y u (only the 
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measured carriers are estimated, as they are the only ones carrying data, thus N v = N F ) 
and, for PSAM, it is H = FF uh y u (the whole channel is estimated). Note that in the 
spectral shaping systems case, again a projection of said captured reference tones on the 
column space of a partial N F by N H FFT matrix is involved. 

The invention can be casted as a Constrained Least Squares estimator by 
expressing the reduced channel length in the following equation : 


Ci 


L yi . 


(10) 


where Cu and CI are parts of an inverse FFT matrix. From this equation, we can write a 
constraint equation as 


(11) 


leading to the Constrained Least Squares solution 


H* = e(C«f^fC f Cf (C f Cf fCuYu 


(12) 


where P s is an orthogonal projection matrix (hermitian and P s 2 = Ps ) which can be 
further simplified to P s = C„ # P C i C u . However, finding an explicit expression for S is 
not straightforward. When N c -N u > N he (otherwise P s is not full rank and the CLS 
solution is not defined) and taking into account that C u F uh + G F| h =0 , the following 
equations show that P s = PFuh . Indeed, the combined projection on S and F uh is shown to 
be a projection on F„h : 
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Hence, the two spaces are the same and the estimators are also the same. As H = F(F uh 
F„h) " ! Fuh y u , the channel estimator is the cascade of a partial IFFT, a weighting matrix 
and a partial FFT. For an arbitrary number of pilots ( > N he ) , the global estimation 
5 scheme (i) to perform a partial initial IFFT, partial as only part of the carriers are 
measured, (ii) weighting the non-trivial part of the channel impulse response width by 
(F U h H Fuh)" 1 and (iii) transforming back to the frequency domain. It is an important aspect 
of the invention to recognize that in the general case of non-equidistant carriers and/or 
when not all pilots or reference tones are present a weighting matrix, being essentially 
1 0 different from the identity matrix is needed in order to preserve optimality . 

Figure 11 shows a flowchart of a three step approach. In an embodiment of the 
invention the step (600) performs Fuh, the step (610) performs (F uh H Fuh)' 1 while step 
S (620) performs F. Said embodiment can be implemented by performing explicit matrix 

K multiplications or for the steps (600), (610) using IFFT and FFT circuitry. 

H= 1 5 Figure 1 3 shows a flowchart of a one step approach wherein reference tones (500) 

"j are transformed into a channel estimate (530) in a single step (900). In one embodiment 

H " said single step approach can execute either Ps from equation (12). In an alternative 

0 embodiment the operation defined in equation (14) is performed or a mathematical 

T, J equivalent, as defined via equation (13). 

O 20 Indeed, it is only if all pilots are present, or if they are regularly spaced, that the 

C method boils down to going from the frequency domain to the time domain, force the 

time channel estimator to be of length N he and going back to the frequency domain as 
shown in [G.G. Raleigh and V.K. Jones. "Multivariate Modulation and Coding for 
Wireless Communication," IEEE Journ. on Special Areas in Communications, 17(5):851- 
25 -860, May 1999.[, wherein no weighting is used. 

The classical ML solution can be applied to a combination of PSAM and 
decision-feedback. Indeed, suppose we use the pilot symbols along with decisions taken 
on the other carriers, then r= h+ v remains valid, with a given Cw , which leads to H = 
P Fhe y . Hence, if the designer can afford the increment in complexity, combination of 
30 PSAM and decision-feedback is desirable. Indeed, figure 7 shows that a difference in 
performance of 2-3 dB can be expected between an all-pilot system (which is equivalent 
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to combined PSAM/DF if decision errors are neglected) and a PSAM system with 8 pilot 
carriers. 

The complexity of the ML estimator is significantly lower than SVD approaches, 
both for spectral shaping and PSAM systems. This low complexity relies on the time- 
5 frequency interpretation and pruning of the (I)FFTs. By construction, Ppuh is a low rank 
matrix (of rank N hee ). Taking its hermiticity into account, it can be written as 

EV =W H (14) 

where V is a matrix of size N c x N he that can be precomputed. Hence, the complexity for 
10 computing the estimator is 2 N„ x N he complex multiplications for the global ML 
estimator. 

Figure 12 shows a flowchart of the two step. In an embodiment said first step 
(700) executes the operation defined by the Hermitian inverse of V while the second step 
(720) performs the operation defined by V itself. 

15 This complexity is about the same as for the SVD approach, however, while [0. 

Edfors, M. Sandell, J.J. van de Beek, S. K. Wilson, and P. O. Borjesson. "OFDM 
Channel Estimation by Singular Value Decomposition," IEEE Trans, on 
Communications}, 46(7):931--939, July 1998.] uses an approximation, (14) is an exact 
maximum likelihood estimate. Further complexity reduction can be obtained by using 

20 the time-frequency interpretation. Indeed, the projection operation can be expressed by 
the cascade of two partial FFTs, weighted by a N he xN he non-identity matrix. With a 
radix-4 implementation of the FFT, the complete estimator would require 1.5 N c 
(Iog 4 (N c )-l) +N he 2 complex multiplications. Furthermore, some additional complexity 
gain can be achieved by using FFT pruning or transform decomposition. Such techniques 

25 lead to a significant gain for the Fourier Transforms. However, the last term ( N he ), due 
to the weighting matrix ( (F^Fun)" 1 ), remains unchanged (figures 3,4). 

When using Pilot Symbol Assisted Modulation, a comb spectrum (figure 2) has to 
be measured, and only the teeth of this comb are used for the FFTs. This particular case 
has been studied by [S. He and M. Torkelson, Computing Partial DFT for Comb 

30 Spectrum Evaluation, IEEE, 1996, Vol. 3, no. 6, June, p 173-175]. In this case, the DFT 
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can be computed with N c /4 +N he /2 +Iog 2 N he - N he complex multiplications, which 
represents a large gain for a large number of carriers. Figures 3 and 4 show complexity 
evaluations of the different approaches: the SVD-based approach of [O. Edfors, M. 
Sandell, J.J. van de Beek, S. K. Wilson, and P. O. Borjesson. "OFDM Channel 
Estimation by Singular Value Decomposition," IEEE Trans, on Communications, 
46(7):93 1-939, July 1998.] and the frequency-time approaches according to the 
invention with plain FFTs, with FFT pruning and with FFT optimized for a comb 
spectrum. 

The complexity for FFT-based solutions is much lower than for the SVD-based 
approach, both for spectral shaping and PSAM systems. Furthermore, simulations show 
that the ML algorithm can work with a significantly smaller N he than the LMMSE, which 
results in a still larger gain than appears in figure 3. For a relatively large number of pilot 
carriers, the main contribution to the complexity is due to the weighting matrix (see the 
N he 2 curves in figures 3 and 4). However, for pure PSAM with regularly spaced pilot 
carriers, it can easily be shown that the weighting matrix ( (Fui^Fuh)" 1 ) is proportional to 
the identity matrix, and complexity is even lower. This special case of our algorithm is 
the frequency correlation part of the algorithm developed by Raleigh and Jones [G.G. 
Raleigh and V.K. Jones. "Multivariate Modulation and Coding for Wireless 
Communication," em IEEE Journ. on Special Areas in Communications}, 17(5):851- 
860, May 1999.], which require regularly spaced pilot carriers or reference tones, while 
the invention is not limited thereto. At the contrary the invention is specially constructed 
for modern telecommunication standards with non-regularly spaced carriers. 

To evaluate the performances of the ML estimator, and compare it with the 
LMMSE algorithm, we simulate a spectral shaping system and a PSAM-based system in 
an indoor radio channel. Two OFDM schemes with 64 and 256 carriers are considered, 
both with QPSK modulated carriers and a cyclic prefix of 16. The 64-carriers scheme is 
simulated with the Hiperlan II zero carriers and the 256-carriers with PSAM. The data 
rate is 25 Msamples/second over the air (i.e. including the cyclic prefix) with a carrier 
frequency of 5.6 GHz. We consider a collection of 120 indoor office-like channels. The 
channel is modeled by means of a ray-tracing technique, considering 20 emitter locations 
and 6 receiver locations, as sketched in figure 5. From the ray-tracing results, it appears 
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that the channel length is of the order of 4 to 6 and can be modeled as having an 
exponentially decaying power delay profile with normalized time constant 2 (for 
simulation of [Edfors et al.]). For the spectral shaping system, the Bit Error Rate (BER) 
is simulated for both LMMSE [O. Edfors, M. Sandell, J.J. van de Beek, S. K. Wilson, and 
P. O. Borjesson. "OFDM Channel Estimation by Singular Value Decomposition," IEEE 
Trans, on Communications}, 46(7): 931 --93 9, July 1998.] and ML estimators and for N he 
ranging from 4 to 16. BER based on exact channel knowledge and raw measurements are 
evaluated for comparisons. Simulation results (figure 6) clearly show that the LMMSE 
suffers from a threshold effect at high SNR. To obtain similar performances for both 
algorithms, N he must be 2 to 4 times larger for LMMSE than for ML. The behavior of 
the ML (and to some extent of the LMMSE) estimator at low SNR gives some insight on 
the influence of N he . Going from higher to lower N he leads to better performance (figure 
6) ; indeed, the values of the channel impulse response beyond the lower N he are below 
the noise level, so that their estimation introduces more noise than relevant information 
about the channel. For the PSAM-based system, we have evaluated both estimators with 
8 and 32 tones and N he = 8 (figure 7). In this simulation, "raw measurements" mean 
measurements on all pilots. As for the spectral shaping system, the flooring effect of the 
LMMSE estimator is essentially limiting it's effectiveness at high SNR's. Noteworthy, the 
ML estimator based on 8 over 256 tones gives similar performance as raw measurements 
on all carriers. 

A low complexity Maximum Likelihood (ML) OFDM channel estimator is 
proposed. It relies on a deterministic model (i.e. no statistical information on the channel) 
that takes the finite delay spread of the channel into account, which is linked to the 
frequency correlation of the channel. Our ML estimator can be interpreted as a translation 
of some initial estimate of the frequency response of the channel to the time domain, 
followed by a linear transformation of this channel impulse response, and retranslation to 
the frequency domain. This interpretation leads to low complexity algorithms, derived by 
combining the partial (I)FFTs involved and a small weighting matrix. Comb spectrum 
due to Pulse Symbol Assisted Modulation allows further pruning of the FFT's. Based on 
theoretical grounds that shade a new light on former solutions, the proposed algorithms 
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have a significantly lower complexity than the low-rank approximation of the LMMSE 
estimator, while being optimal. 

The concept of channel estimation proposed here fits within telecommunication 
devices as described in US 09505228. Figure 8 shows an architecture of a signal 
5 receiving circuit, wherein said channel estimation method fits, a fast Fourier transform 
circuit (2), a symbol reordering block (3), a time synchronization circuit (4), performing 
carrier and clock frequency offset estimation, an equalization circuit (5), a demapping 
circuit (6) and a channel estimation circuit (7), executing the described channel 
estimation method. Said channel estimation circuit provides channel information (8) for 
10 said equalization circuit (5). Said circuit can easily be used also for signal transmitting 
when a mapper circuit (1) is added. 

In an embodiment said time synchronization circuit (4) determines control 

BSW. 

3 information on that part of its input, being a non-orthogonal frequency division 

H multiplexing signal, for instance partly comprising of a Pseudo-Noise sequence. Said 

=* 15 control information can be signal level information, carrier offset or timing information. 
lj The signal conversions in the time synchronization circuit (4) comprise at least of carrier 

offset compensation and guard interval removal. 
3 In an embodiment said mapper circuit (4) is a serial to parallel converting circuit. 

*j In an embodiment said symbol reordering circuit (3) performs guard interval 

3 20 introduction and symbol reordering for the modulating dataflow and symbol reordering, a 
2 first step for despreading and zero carrier grouping for the demodulating dataflow. 

In an embodiment said equalization circuit (5) performs corrections for channel 
distortions, thereby using the estimated channel information (8), on signals, said 
corrections are preferably complex multiplications or divisions. 
25 In an embodiment said demapping circuit (6) extracts bits from the carriers of the 

frequency domain representation of the equalized signals. 

The inputs to said receiver typically go though the time synchronization circuit 
(4), then though the fast Fourier transform circuit (2) and thereafter through the symbol 
reordering circuit (3). The symbol reordered signal is then fed into the equalization circuit 
30 (5), whereafter the demapping circuit (6) is used for extracting the final output. The 
mapper circuit (1) when available is connected to the fast Fourier transform circuit (2). 
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It is important to note that channel estimation must not be confused with carrier 
and clock offset estimation. More in particular, in the telecom device shown, said 
determining at and in the receiving peer happens on said reference tones, being corrected 
for carrier and clock frequency offsets. Hence said channel estimation happens after 
5 carrier and clock frequency offset correction. Note that Figure 1 shows a representation 
of an OFDM system, wherein the telecom devices, shown in Figure 8, fit. Indeed said 
devices can be used both at the transmitting (left side of Figure 1) and the receiving side 
(right side of Figure 1). The information (8) provided by the channel estimation circuit 
(7) are the channel value h[0], ...h[N-l], with N here the amount of carriers in the 
1 0 channel. The equalization boils down to the complex divisions, as shown in Figure 1 . The 
concept of channel estimation proposed here uses frequency pilots to estimate the 
frequency response times and is based on a particular non-iterative linear estimation of 
5 the channel responses, further enabling frequency domain channel equalization. The 

?? derived least mean square approach explicitly exploits the fmiteness of the channel 

15 impulse response. 

y In an aspect of the invention essentially digital circuits or devices for determining 

=b a maximum likelihood frequency domain estimation of the channel response of a channel 

3 between at least one transmitting peer and at least one receiving peer, according to the 

jj invented method, are presented. 

2 20 In a first embodiment a device at a receiving peer comprising means for capturing 

^ Nu reference tones, transmitted from a transmitting peer to said receiving peer and 

means for determining at said receiving peer from said Nu reference tones said 
maximum likelihood frequency domain estimate of said channel response at N F 
predetermined frequencies by directly exploiting the fmiteness of the time response of 
25 said channel. 

Said means for capturing comprises of circuitry typically found in a receiver, such 
as an antenna, a low noise power amplifier, a bandpass filter, an analog-digital converter. 
In a further embodiment said means include a time synchronization circuit, a FFTYIFFT 
circuit and a symbol reordering block. 
30 Said means for determining said maximum likelihood frequency domain estimate 

comprises of a digital circuit capable of executing the necessary computations and 
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storage devices for storing predetermined and intermediate computated values. Said 
digital circuit can be a general purpose programmable circuit or a dedicated custom 
designed hardware block or a combination of both. 

Figure 9 shows a digital circuit (10) for determining at said receiving peer from 

5 Nu received reference tones (100) said maximum likelihood frequency domain estimate 
(200) of said channel response at N F predetermined frequencies by directly exploiting the 
finiteness of the time response of said channel. Said digital circuit (10) comprises a 
frequency-to-time domain circuit (20) for transforming said received reference tones in 
time domain measurements (110); an estimation circuit (30) for determining from said 

10 time domain measurements maximum likelihood time domain estimates (120); and a 
time-to-frequency domain circuit (40) for transforming said maximum likelihood time 
domain estimates in frequency domain for obtaining said maximum likelihood frequency 
domain estimate. Said estimation circuit (30) comprises a storage device (50), being 
capable of storing a N H by N H weighting matrix, being essentially different from the 

15 identity matrix, with N H the length of said finite time response; and circuitry (60) such as 
a plurality of multipliers for making linear combinations of said time domain 
measurements thereby using coefficients (300) from said weighting matrix. In an 
embodiment said circuits (20) and (40) comprise IFFT and FFT circuits. In another 
embodiment said circuits (20) and (40) also comprise of storage units (70), (80) and 

20 computation circuitry such as a plurality of multipliers (65), (75), using stored 
information (310), (320). Note that said storage devices can define one big memory (90). 
Also part of said multipliers (65), (60), (75) can be shared (91) for executing the 
computations. Said frequency-to-time domain circuits (20) outputs are connected 
logically to said estimation circuits (30) inputs and said estimation circuit (30) outputs are 

25 connected logically to said time-to frequency domain circuits (40) inputs. 

From the formula H = F(F uh H F uh ) A F uh y u , one can derive the operations of the 
frequency-to-time circuit (20), performing Fuh, either via a partial FFT or via 
multiplication with the appropriate matrix, of the estimation circuit (30), performing the 
inverse matrix, via multiplication and of the time-to-frequency circuit (40), performing F. 

30 In said memory (70) in Figure 9, at least the matrix Fuh is stored. In said memory 

(80) in Figure 9, at least the matrix F is stored. In memory (50) at least the non-identity 
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matrix, being the inverse of the product of the Hermitian inverse of Fun and Fuh, is 
stored. Potentially said matrices are stored in a single memory structure (90). Said banks 
of multipliers (65), (60), (65) can be a single bank (91) of multipliers or at least share 
some of the multipliers. 

5 Figure 10 shows a digital circuit implementing a two-step approach comprises a 

first subcircuit (1020) transforming said received reference tones (100) into an 
intermediate results (400); and a second subcircuit (1040) transforming said intermediate 
result (400) into said maximum likelihood frequency domain estimates (200). Said first 
subcircuits (1020) outputs, carrying said intermediate results (400) are logically 

1 0 connected to said second subcircuits (1 040) inputs. 

From equation (14) one can derive the operation of the first subcircuit (1020), 
performing the operation of Hermitian inverse of V and the second subcircuit (1040), 
performing the operation of V itself. 

In an embodiment said circuit for transforming said reference tones in time 

1 5 domain measurements comprises a plurality of multipliers for multiplying a partial N H by 
Nu IFFT matrix with said reference tones with N H the length of said finite time response. 

In an embodiment said circuit for transforming said reference tones in time 
domain measurements comprises a partial N H by Nu IFFT transform circuit with N H the 
length of said finite time response. 

20 In an embodiment said circuit for transforming said maximum likelihood time 

domain estimates into said maximum likelihood frequency domain estimate comprises a 
plurality of multipliers for multiplying a partial N F by N H FFT matrix with said maximum 
likelihood time domain estimates. 

In an embodiment said circuit for transforming said maximum likelihood time 

25 domain estimates into said maximum likelihood frequency domain estimate comprises a 
partial N F by N H FFT transform circuit. 

In an embodiment said circuit for determining from said time domain 
measurements maximum likelihood time domain estimates comprises a storage device, 
being capable of storing a N H by N H weighting matrix, being essentially different from 

30 the identity matrix, with N H the length of said finite time response and circuitry for 
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making linear combinations of said time domain measurements thereby using coefficients 

from said weighting matrix. 

In an embodiment said digital circuit (10) for determining at said receiving peer 
from Nu received reference tones said maximum likelihood frequency domain estimate of 
said channel response at N F predetermined frequencies comprises a plurality of 
multipliers for multiplying a N F by Nu matrix with said reference tones. 

Note that power-optimized FFT circuitry, as disclosed in US 09613181, can be 
used for said FFT and IFFT circuits needed for the channel estimation and preferably also 
for the FFT block in the reference tone capturing means. 
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